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The final parsec problem: aligning a binary with an external 
accretion disc 
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ABSTRACT 

We consider the interaction between a binary system (e.g. two supermassive black holes or 
two stars) and an external accretion disc with misaligned angular momentum. This situation 
occurs in galaxy merger events involving supermassive black holes, and in the formation of 
stellar-mass binaries in star clusters. We work out the gravitational torque between the binary 
and disc, and show that their angular momenta Jh,Jd stably counteralign if their initial orienta- 
tion is sufficiently retrograde, specifically if the angle 9 between them obeys cos 6 < -/(j/2/b, 
on a time short compared with the mass gain time of the central accretor(s). The magnitude 
7b remains unchanged in this process. Counteralignment can promote the rapid merger of su- 
permassive black hole binaries, and possibly the formation of coplanar but retrograde planets 
around stars in binary systems. 
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1 INTRODUCTION 

Galaxy mergers are commonly thought to be the main mechanism 
driving the coevolution of galaxies and their central supermassive 
black holes (SMBH). In such a merger we expect the forma- 
tion of a SMBH binary in the centre of the merged galaxy. 
Gravitational waves quickly drive the binary to coalesce if the 
orbital separation can be shortened to < lO^^pc. The binary 
may stall at a separation greater than this if the interaction with 
the merged galaxy is not efficient enough in extracting orbital 
angular momentum and energy. For the stellar component of the 
galaxies this occ urs at approximately a parsec, creating "the final 
parsec problem" l IMilosavlievic & Merritt. 20011) . There have been 
many papers exploring potential solutions to this problem, for 
example a sling-shot mechanism involving a triple SMBH system 
jlwasawa et al.. 200q) . efliicient refilling of the binary loss cone by 
angular momentu m exchange between stellar orbits and a triaxial 
dark matter halo l IBerczik et al.. 200q ) and also the evolution of 
the binary with a prograde a ccretion disc (circumbinary discs : 
Armitage & Nataraian 20051; MacFadven & Milosavlievic 20081 ; 



Lodato et al. 20091 ; Cuadra et al. 2009| and emb edded discs: 
Escala et al. 20051; IPotti et al. 200iDotti et al. 20091) . In a recent 



paper l INixon et al.. 20 111) we explored the evolution of a binary 
interacting with a retrograde circumbinary accretion disc. We 
showed that this is more efliicient than a prograde disc in removing 
binary orbital angular momentum and energy. This is simply 



because there are no orbital resonances between the binary and the 
disc and thus there is direct accretion of retrograde gas onto the 
binary. 

Here we consider the alignment process between a binary sys- 
tem and an external misaligned accretion disc. This situation can 
arise in at least two astronomical contexts. First, a merger event 
between galaxies can produce a SMBH binary in the centre of the 
merged galaxy, and this or a later accretion event may surround the 
hole with a disc of accreting gas. A similar situation arises during 
the formation of stars in a cluster A binary system may form, but 
also capture gas into an external disc. 

In both of these cases, there is no compelling reason to as- 
sume that the binary and disc rotation are in itially parallel or even 
roughly coaligned (cf iKing & Pringle 200q) . As we shall see, the 
gravitational interaction between the binary and the disc generates 
differential precession in the disc gas, and thus viscous dissipation. 
This gives a dissipative torque which vanishes only when the binary 
and disc angular momenta /b,/d are either parallel or antiparallel. 
In all s uch cases, the torque diffuses the tilt or warp thro ugh the 
disc (cf lPringle 199 JPringle 199^ ; lwiiers & Pringle 19991) driving 
the system to one of these equilibria. The existence of a warp makes 
the precise definition of disc angular momentum /j quite subtle and 
we return to this point in the Discussion. 

The binary-external disc interaction is very similar to the 
effect of the Lense-Thimng (LT) p recession on an accretion 
disc around a spinning black hole JBardeen & Petterson 19751 ; 
Pringle 19921; [Sch t 
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iNelson & Papaloizou 200ol : iLodato & Pringle 20061 etc) if we 
replace J^, by the hole spin angular momentum /i,. For some 
years it was thought that the LT interaction always led to co- 
alignment (i.e. /i, and J^ parallel). However iKing et al. (20051) 
(hereafter KLOP) showed on very general grounds that coun- 
teralignment does occur, if (and only if) the initial angle 6 
between J^ and Jh satisfies cos 8 < - ](^j2]^, where J^ = \Ji\ 
and A = L^hl- IScheuer & Feller (199q) had implicitly assumed 
7d » ^h a nd so enforced co-alignment. With this restric- 
tion lifted, iKing & Pringle (200dKing & Pringle (2007h and 



iKing et al. (20081) showed that accretion from a succession of 
randomly-oriented discs leads to spindown of the supermassive 
black hole, allowing rapid mass growth. 

In this paper we examine the alignment process for a bi- 
nary and an external disc. We show that the argument of KLOP 
is generic, and that the disc and binary counteralign if and only 
if cos 9 < -Jil2J\,. As a result it is quite possible for SMBH 
binaries to be surrounded by a co mpletely retrograd e disc which 
strongly promotes coalescence (cf iNixon et al. 20111) . In the case 
of a newly-formed stellar binary, the presence of a counteraligned 
disc can lead to the formation of planets with retrograde orbits. 



2 THE BINARY-DISC TORQUE 

We consider a binary system with masses Mi , Mj and a circular 
orbit, with the binary angular momentum vector pointing along the 
z-axis of cylindrical polar coordinates (R, i/>, z). For simplicity we 
assume M2 « Mj and place Mi at the origin, with M2 orbiting at 
radius a in the {R, ^) plane (our conclusions are not affected by this 
assumption). The orbit has angular velocity 



nb = 



G(M| + M2) 



(1) 



Now we consider a disc particle in an orbit about the binary 
at radius /? » a. If both the small quantities M2/M1 and alR ac- 
tually vanished, the particle's orbit would be a circle, with angu- 
lar velocity (GMi/i?)''^. When these quantities are small but finite 
they induce various perturbations in the orbit. Some of these per- 
turbations have (inertial-frame) frequency IQ.^, and higher multi- 
ples. These are oscillatory, and have no long-term secular effect. 
Long-term effects on the orbit, and hence eventually on the disc, 
come from the zero-frequency (azimuthally symmetric m = 0) 
term i n the binary pote ntial. This point is considered in more de- 
tail in lSate et al. (20001) , who considered the related problem of a 
disc around the primary mass Mi (i.e. Mj ^ Mj , but R <g: a). 

Physically this in = term is given by replacing the orbiting 
mass M2 with the same mass spread uniformly over its orbit, i.e. 
a ring of mass M2 and radius a in the (R, if}) plane. Adding in the 
potential from the fixed point mass Mi at the origin we find the 
effective gravitational potential experienced by a disc particle as 



<S>(R,z) = -- 



GMi 



GM2 r'" dcp 

2k Jo r 



(R^ + z^y- 2n ' ^^^ 

where r is the distance between the particle position and a point on 
the ring at (a, (p, z), i.e. 



■ R' + a- + z - 2Ra cos c 



(3) 



We now expand equation[2]in powers ofa/R and z/R, keeping terms 
only up to second order This gives 



<i>{R,z) 



G(M, + M2) GM2a^ G(M, + Mj)^^ 



9GM2fl^z^ 

8^5 



(4) 



The orbital frequency Q of the particle subject to this potential is 
given by 



RdR 

and its vertical oscillation frequency v by 

dz'' 

both evaluated at z = 0. The nodal precession frequency is Q.p 
Q.- V and we find 

nl/2 



y2 



(5) 



(6) 



njR) = - 



G(Mi + M2) 



R^ 



M, 



Ml -H M2 /?2 ■ 



(7) 



This frequency is very similar to that for LT pr ecession around 
a spinning black hole (e.g. IScheuer & Feller 1996h . which goes as 
R^^ rather than the i?"^'^ here. Equation [7] i s formally almost iden- 
tical to the precession frequency found bv lSate et al. (2000l) for a 
disc around the primary, although derived for /? » a rather than 
R ^ a. The same argument as in that paper shows that if the 
disc and binary axis are misaligned by an angle 6 (called S in 
ISate et al. 20000 with < 6 < 7t/2, the precession frequency is 
just multiplied by cos 6. The opposite case with the disc somewhat 
counteraligned (i.e. 9 > n/2) is equivalent to the 6 < n/2 case with 
the binary angular momentum reversed. But this reversal leaves the 
precession frequency unchanged, since we are dealing only with 
the m = part of the potential. So for all 9 with < 9 < n the 
precession frequency is 



QJO) = Qr, |cos e\ . 



(8) 



R 



4R3 



2R3 



This result differs from the LT case, where the factor cos 6 appears 
without modulus signs. 



3 CO- OR COUNTER-ALIGNMENT? 

We have shown above that the effect of the binary potential on the 
disc is to induce precession of the disc orbits. This precession is 
strongly dependent on radius: rings of gas closer to the binary pro- 
cess faster. The differential precession creates a dissipative torque 
between adjacent rings of gas tending to make 9 — > 0, ;r so that 
the precession ultimately vanishes. 

The precession timescale in the disc increases with radius 
(cf. equation ITJI. The torque therefore acts faster at smaller radii 
to CO- or counteralign disc orbits with the binary plane. This 
leads to the creation of a warp in the disc, where the inner parts 
are co- or counteraligned and the outer parts are still misaligned 
(cf. Fig. [TJ. This warp propagates outwards, eventually co- or 
counter aligning the entire disc with the binary plane. This effect 
was solved numerically for discs warped under the LT effect by 
ILodato & Pringle (2006). 

Now we argue as in KLOP that since each ring feels a preces- 
sion, the resultant back-reaction on the binary is a sum of preces- 
sions, which is just a precessio n. This argument is equivalent to that 
presented in lBate et al. (2000l) who argue that because the binary 
potential is symmetric about the plane of the binary, the disc-binary 
torque cannot have a component in the direction of /b. Accordingly 
Ji, can only precess. These arguments show that we can write the 
torque on the binary in the same form as the LT-induced torque on 
a spinning black hole considered in KLOP, i.e. 
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^ = -KAh A/d] - K2[h A (/b A/d)]. (9) 

Here K[, K2 are coefficients depending on disc properties. The first 
term gives the magnitude and sign of the torque inducing the pre- 
cession. It does not change the aUgnment angle 0. The second term 
describes the torque which changes 9. The same arguments as in 
KLOP for the LT case, and lSate et al. (20001) for a disc around the 
primary, show that dissipation in the disc requires Ki to be a posi- 
tive quantity. Its magnitude depends on the properties of the disc 
and the binary. The one difference from the LT case is that the 
|cos d\ dependence means that the sign of the coefficient A", can 
be either positive or negative. But this difference has no effect on 
the conditions under which the disc and binary co- or counteralign. 
These are formally identical with the ones for the LT case derived 
by KLOP, with the binary angular momentum /b replacing the hole 
spin angular momentum /i,- The process obviously has a different 
timescale specified by the different magnitude of K2. 

The same arguments as in KLOP now show that the magnitude 
^b of the binary angular momentum remains constant, while the di- 
rection of /b aligns with the total angular momentum J, = Jb +Ji, 
which is of course a constant vector. During this process the magni- 
tude of jj decreases because of dissipation (KLOP). Counteralign- 
ment (0 — > k) occurs if and only if J^ > Jf. By the cosine theorem 



Jbi 



jj = JI + J^- 2/b-/d cos(;r - 0), 
so this is equivalent to 

COS 6 < . 

2yb 



(10) 



(11) 



Thus counteralignment of a binary and an external disc is possible, 
and requires 

> n/2, Ji < 2Jb. (12) 



4 DISCUSSION 

So far in this paper we have avoided fully spelling out the mean- 
ing of the disc angular momentum /d- This is complicated because 
the binary torque falls off very strongly with radius, and so a large 
contribution to the angular momentum in a distant part of the disc 
may be irrelevant to the al ignment process, or affe ct this process in 
a time-dependent way (cf lLodato & Pringle 20061) . Sections 3 and 
4 of KLOP discuss these questions in more detail. Effectively /d 
can be thought of as the disc angular momentum inside the warp 
radius, and therefore a time-dependent quantity. 

At early times J^ is small, as only a fraction of the total gas 
interacts with the binary. Counter-alignment may occur if (? > ;r/2, 
but at later times, as yd grows and more gas is able to interact 
with the binary, alignment eventually happens (when 7d > 2yb)- 
So if S > n/2, even for yd > 2yb we expect ~ 2yb|cos6| of 
disc angular momentum to counteralign with the binary before 
the outer disc comes t o dominate and enforce coalignment (cf. 
iLodato & Pringle 20061) . 

The typical timescale for co- or counter-alignment for a 
SMBH binary is 



^bini 






(13) 



where i?„ is the warp radius, Jd(Ry,) is the disc angular momen- 
tum within if„, and V2 is the vertical disc viscosity. This is identi- 
cal to the formal expression for LT alignment of a spinning black 
hole if we replace the spin angular momentum J^ with yb (cf 



CO- or counter-aligned 
inner disc 




warp region 



misaligned outer disc 

Figure 1. The warped disc shape expected after the inner disc co- or 
counter-aligns with the binary plane but the outer disc stays misaligned. 
Eventually the entire disc will co-or counter-align with the binary plane, 
depending on the global criterion (egn lllt . Note that in practice precession 
makes the warp non-axisymmetric. 



IScheuer & Feiler 19960 . The warp radius is given by equating the 
precession time l/QpiR) to the vertical viscous time R^/v2 . Inside 
this radius the precession timescale is short and the disc dissipates 
and CO- or counter-aligns with the binary plane. Outside this ra- 
dius the disc is not dominated by the precession and so maintains 
its misaligned plane. The connecting region therefore takes on a 
warped shape shown in Fig.[T] As time passes the warp propagates 
outwards and co- or counter-aligns the entire disc with the binary 
plane. 

Approximating the disc angular momentum as 



J,(R„) ~ ttRIKGMR^) 



1/2 



(14) 



with E the disc surface density and M = Mi + M2, and using the 
steady-state disc relation M = 3nvll we find 



^bini 



V2 M' 



M,\R. 



where we have also used 



, Ga 

A = MiM2\ — 

M 



1/2 



(15) 



(16) 



Since vi < V2 jPapaloizou & Pringle. 19831) . a ^ R„ and M2 < Mi, 
we see that alignment takes place on a timescale shorter than the 
mass growth of the central accretor(s). 

The timescale il5b is directly analogous to the expression 






(17) 



for alignment under the LT precession, where a, < 1 is the Kerr 
spin parameter and R^ the Schwarzschild radius of the spinning 
hole. Evaluating R„ in the two cases we find 

(18) 



tu 



fbir 



31^ 

2 



[M2/M1J [rj 



Thus in general, provided we assume that the ratio vi/v2 is similar 
in the two cases and that the hole spin is not rather small (a, < 
(i?s/a)''^(M2/Mi)), then the binary-disc alignment is rather faster 
than the corresponding process for spinning black holes. 

Our result has significant consequences for SMBH bi- 
naries. For random orientations, equation [12] shows that ini- 
tial disc angles leading to alignment occur significantly more 
frequently than those giving counteralignment only if yd > 
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2Jf,. (In the LT case this fact leads to a slow spindown 
of the hole, becau se retrograde ac cretion has a larger ef- 
fect on the spin. iKin g et al. 2008.) A number of studies 
I Armitage & Nataraian 2005: MacFadven & Milosavlievic 20081 : 



ICuadra et al. 2003 : iLodato et al. 20090 have shown that prograde 
external discs are rather inefficient in shrinking SMBH binaries and 
solving the last parsec problem. This is essentially because of res- 
onances within the disc. In contrast, the slightly rarer retrograde 
events have a much stronger effect on the binary. These rapidly pro- 
duce a counterrotating but coplanar accretion disc external to the 
binary, which has no resonances. We note that lNixon et al. (201 ih 
show that the binary gradually increases its eccentricity as it cap- 
tures negative angular momentum from the disc, ultimately coa- 
lescing once this cancels its own. A non-zero binary eccentricity 
changes the detailed form of the perturbing potential from that in 
equation|4l but cannot change the precessional character leading to 
the torque equation ^. Our results remain unchanged, particularly 
the counteralignment condition l ll2t . apart from minor modifica- 
tions of the timescale l ll5t . 

Thus in a random sequence of accretion events producing ex- 
ternal discs, the prograde events have little effect, and the retrograde 
ones shrink the binary. In particular, a sequence of minor retrograde 
events with Ji < Ji, has a cumulative effect and must ultimately 
cause the binary to coalesce once the total retrograde 2 ■^d = -^b- 
This is important, since the disc mass is limited by the onset of self- 
gravity to Md < (H/R) Ml (cf King, Pringle and Hofmann 2008). 
Coalescence will then occur once the retrograde discs have brought 
in a total mass M2, i.e. once a sequence of > {M2/Mi){R/H) retro- 
grade discs have accreted. For minor mergers this requires at most 
a few randomly oriented accretion disc events, rising to a few hun- 
dred for major mergers (q > 0.1). 

We note finally that similar considerations apply in planet- 
forming discs aro und stellar binary systems, which can also be ini- 
tially misaligned JBate et al.. 20100 . This may offer a way of mak- 
ing re trograde planets in bina ries, as recently suggested for v Oc- 
tantis JEberle & Cuntz. 2010h . 
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